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ABSTRACT 

Recent observations reveal that spectral breaks at ^GeV are commonly present in 
Galactic 7-ray supernova remnants (SNRs) interacting with molecular clouds and that 
most of them have a spectral (E 2 dF/dE) "platform" extended from the break to lower 
energies. In paper I (Li & Chen 2010), we developed an accumulative diffusion model 
by considering an accumulation of the diffusive protons escaping from the shock front 
throughout the history of the SNR expansion. In this paper, we improve the model 
by incorporating finite- volume of MCs, demonstrate the model dependence on particle 
diffusion parameters and cloud size, and apply it to nine interacting SNRs (W28, W41, 
W44, W49B, W51C, Cygnus Loop, IC443, CTB 37A, and G349.7+0.2). This refined 
model naturally explains the GeV spectral breaks and, especially, the "platform" s, 
together with available TeV data. We find that the index of the diffusion coefficient 
8 is in the range of 0.5-0.7, similar to the galactic averaged value, and the diffusion 
coefficient for cosmic rays around the SNRs is essentially two orders of magnitude lower 
than the Galactic average (\ ~0.01), which is a good indication for the suppression of 
cosmic ray diffusion near SNRs. 
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1 INTRODUCTION 

Supernova remnants (SNRs) are commonly believed to be 
one of the most important acceleration sites for the cosmic 
rays (CRs) below the "knee" in our Galaxy. Diffusive shock 
acceleration (DSA) is the prevailing acceleration mechanism, 
which can naturally gain a power-law population of relativis- 
tic electrons/protons (Blandford & Eichler 1987). Although 
multi-wavelength analysis has been intensively made for the 
emission of SNRs, it is still in hot debate whether the 7- 
rays from them are of hadronic or leptonic origin. Recently, 
7-ray telescopes of new generation, such as Fermi-LAT and 
H.E.S.S., provided more and more clues that SNRs inter- 
acting with molecular clouds (MCs) may emit 7-rays arisen 
from 7T° decay via proton-proton collision (e.g., Abdo et al. 
2009, 2010a,b,c,d; Aharonian et al. 2008a,b). 

Interacting SNRs, distinguished by several kinds of evi- 
dence such as OH masers, molecular line broadening, etc. (see 
Jiang et al. 2009 and references therein), represent a promis- 
ing class of 7-ray sources. It is seen that some of the SNRs ex- 
hibit a spectral break at ~GeV, such as W28, W44, W51C, 
IC443, etc. Most of them have a spectral (E 2 dF/dE) "plat- 
form" extended from the break to low energies. By coin- 
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cidence, the 7-ray emitting SNRs which harbor OH masers 
(Hewitt et al. 2009) all have GeV breaks. The spectral break 
can be explained in the hadronic scenario either by different 
acceleration effects in the SNR-MC system, such as Alfvcn 
wave evanescence in the weakly ionized dense gas (Malkov, 
Diamond & Sagdeev 2010), reacceleration model for crushed 
clouds (Uchiyama et al. 2010), and two-step acceleration 
model in the reflected shocks (Inoue, Yamazaki & Inutsuka 

2010) , time-dependent two-zone model (Tang et al. 2011), or 
by diffusion effects of CRs escaping from SNR shock (Li & 
Chen 2010, hereafter Paper I; Ohira, Murase, & Yamazaki 

2011) . 



In Paper I, we developed an accumulative diffusion 
model by considering the accumulation of the diffusive pro- 
tons escaping from the shock front throughout the history of 
the SNR expansion; the power-law distribution is assumed 
for the escaping protons and the spectrum of the diffusive 
protons at any point near the SNR is obtained. This model 
is used to explain the GeV break and GeV/TeV flux dis- 
crepancy for the four sources around SNR W28. Comparing 
with the point source injection (Aharonian & Atoyan 1996; 
Gabici et al. 2009), our model is more physical when ap- 
plied to the finite-size accelerator, such as SNRs, and small 
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distance between the accelerator and MC is allowed^ Al- 
most meantime, Ohira et al. (2011) have also developed a 
diffusion model by considering finite-size MCs interacted by 
the CRs escaping from the SNRs which are imbedded in 
MCs and used the model to fit the 7-ray spectra of four 
SNRs (W28, W44, W51C, and IC443). In their model, the 
time evolution of maximum particle energy is considered and 
mono-momentum is assumed for the escaping particles and 
the spectral shape of the diffusive protons is analytically 
derived. 

In this paper, we improve the accumulative diffusion 
model by incorporating finite-volume of MCs (Sec. [2J and 
apply it to a series of interacting SNRs to reproduce the 
GeV-TeV spectra as well as the "platform" below the break 
at around GeV (Sec. [3]). 



2 ACCUMULATIVE DIFFUSION MODEL 
WITH FINITE- VOLUME MCS 

2.1 Model description 

In Paper I, the distribution function, at an arbitrary point 
(at radius R c ), of the energetic protons that escape from the 
spherical shock front (at radius R s ) throughout the history 
of the SNR expansion is given by 

/cum(£p, Rc, tagc) = / W /; f(E p , RiRcti, 9, <t>),t dil ) 

R%(ti) sin 6 dOdijidti, (1) 

where i age is the age of the SNR, ti the time at which the 
proton escapes from SNR surface, idif the diffusion time after 
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tn 
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distribution function, at point C, of the energetic protons 
that escape from unit area at an arbitrary source point S 
on the spherical shock front surface (see Paper I for relevant 
notations). The diffusion coefficient is assumed to be in the 
form of D(Ep) = 10 28 x(S p /10GeV)' 5 cm 2 s- 1 , where X is the 
correction factor of slow diffusion around the SNR (Fujita 
et al. 2009) and 8 is the energy dependent index of diffusion 
coefficient. The distance between the SNR and MC does not 
need to be larger than the size of SNR, namely condition 
R c 3> R s is not required as opposed to the point source 
injection models. 

Following the treatment of the dynamical evolution 
of SNR in Paper I, we use the Sedov- Taylor law R s = 
(2.026-E SN Rt 2 /p C i) 1/ ' 5 for the adiabatic phase, where p = 
1.4m H w is the density of the interstellar (intercloud) 
medium, and R s = (147ei5s N R_R 2 i 2 /47rp ) 1/ ' 7 for the radia- 
tive phase, where R t is the transition radius from the Sedov 
phase to the radiative phase, -Esnr is the supernova explo- 
sion energy, and e is a factor equal to 0.24 (Blinnikov et al. 
1982; see also Lozinskaya 1992). 

In Paper I, we assume that all the MC mass is concen- 
trated in a point. Here we take the finite-volume of cloud into 
account and in order to describe the problem with spheri- 
cal coordinates, we approximate it as a truncated cone which 



1 Actually, our model in Paper I has been tested for the case 
that the SNR has a very small radius (e.g., for a very high ambi- 
ent density), which can be approximated as a point-like source, 
and reproduced the results of the classical point-like continuous 
injection case using the same model parameters. 



subtends solid angle Q, at the SNR center and has a thickness 
AR C . By averaging the distribution function / cum (Eq. {TJ) 
over the finite-volume, we obtain the mean distribution of 
the energetic protons in the MC at f ago as a function of E p 
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' R c -AR c /2 

In particular, for Q — 4tv, we have the case that the SNR is 
enclosed by a molecular shell (as also assumed in Ohira et 
al. 2011). 

In the calculation of the 7-rays from the nearby MC (of 
mass Mc) due to p-p interaction, we use the analytic photon 
emissivity dN^/dE-y developed by Kelner et al. (2006; also 
see Paper I). The gained photon energy is generally 10% of 
the energetic proton energy (e.g. Katz & Waxman 2008). 



2.2 Model Performance 

We now explore the model parameters and show how the 
spectra of escaping protons depends on the parameters. In 
the model described above, for a certain SNR, the model 
parameters are of two sorts: particle diffusion (p, 8, x) and 
MCs {Rc, AR). 

Figure 1 shows the dependence of the spectral shape 
of the energetic protons in the finite cloud volume on each 
parameter with other parameters fixed. Here we assume a 
middle-aged SNR with _R S nr = 20pc, t = 27000 yr, n = 
lcm -3 , and E S nr = 10 51 erg and a proton illuminated MC 
at R c = 25pc. 

The spectral dependence on the MC thickness AR C is 
presented in Figure la, where we adopt \ = 0.01, p — 2.2, 
and 8 = 0.5. We plot three curves for different thicknesses: 
ARc -¥ 0, AR C = 8 pc, and AR C = 10 pc. In the AR C 
case, the cloud volume is not considered, similar to the treat- 
ment in Paper I. We see that there are spectral breaks in 
the diffusive protons and find that, with the thickness in- 
creased, the slope change around the break becomes gentle 
and a "platform" -like pattern appears and gets broadened. 
This pattern is formed by the combination of particle dis- 
tribution at different diffusion distances (within the cloud 
thickness) , as the spectral peak shifts to higher energies with 
the increase of distance. The AR C = 10 pc case represents 
the situation that the MC is in contact with the SNR shock 
surface (R c — R s + AR c /2; such case is called 'contact' 
case hereafter) , in which the proton spectrum does not show 
a low-energy cutoff/turn-over since some of the low-energy 
protons can diffuse into the cloud. 

The spectral dependence on the escaping protons energy 
index p is given in Figure lb, with parameters x = 0.01, 
AR = 15pc, and 8 — 0.5 used. Three values of p are used for 
comparison: p = 1.8 is the particle spectral index predicted 
by the DSA mechanism considering non-linear effects, p — 2 
is the typical index of Fermi-type accelerated particles, and 
p = 2.2 is the source spectral index of the Galactic CRs 
derived from the diffusion model GALPROP (Strong et al. 
2007). For the 'contact' case (A_R C = 10 pc), it can be seen 
that the slopes (both below and above the break energy) 
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of the particle distribution increase with p. For the 'non- 
contact' case (R c > R s +AR c /2, e.g. AR C = 8 pc), the slopes 
above the break are very similar to those of the 'contact' 
case, and the low energy cut-off becomes more prominent 
with the increase of p. 

Figure lc shows the spectral dependence on the diffu- 
sion coefficient index 5 (with \ = O.Of, AR = f5pc, and 
p — 2.2 fixed). The value of 5 is usually adopted in the 
range of 0.3-0.7 (Berezinskii et al. (f990)). Here we use 
three values of 8: S =0.3, very close to the Kolmogorov type 
diffusion index (0.33); S =0.5, for the Kraichnan type dif- 
fusion; and S =0.7. This parameter strongly influence the 
spectral pattern. The slope of the proton spectrum above 
the break conforms the relation 02 ~ p + S for the contin- 
uous point source injection case (Aharonian 2004). For the 
'contact' case (AR C = fO pc), there exhibits a "platform" 
pattern extended to low energy, while for the 'non-contact' 
case (AR C = 8 pc), the low energy cut-off becomes more 
significant with the increase of S. 

As shown in Figure Id, different correction factor of 
diffusion coefficient \ gives different flux and width of the 
"platform". With the increase of the % value, the number of 
the low energy protons increases and the number of the high 
energy protons decreases, namely, the particle distribution 
gets softened. This is because faster diffusing high energy 
protons diffuse into farther distance and the low energy ones 
become dominant at the concerned distance. Meanwhile, the 
width of "platform" becomes narrower and the slope of the 
proton spectra above the break energy keeps unchanged. 



3 APPLICATION TO INDIVIDUAL SNRS 

Extended GeV emission associated with interacting SNRs 
has recently been revealed by a series of Fermi-LAT obser- 
vations (Abdo et al. 2009, 20f0a,b,c,d; Castro & Slane 20f0), 
some of which are also spatially resolved by the H.E.S.S. TeV 
observations. These SNRs form a special interesting class, 
which may be beneficial to unveil the enigma of the origin 
of the Galactic CRs. One of the most conspicuous features of 
such a class is the commonly present, often gentle, spectral 
break at around f-f OGeV, which is difficult to be reproduced 
by leptons via inverse Compton emission (e.g. Abdo et al. 
20f0a). In addition, the TeV spectra are mostly steep, with 
the photon indices sometimes as high as ~ 3.0, which is also 
hard to be explained by the primary accelerated protons, 
but can be done by the proton diffusion scenario (e.g. Aha- 
ronian & Atoyan f996; Torres et al. 2008; paper I; Ohira et 
al. 20fl). Here we apply our improved accumulative diffu- 
sion model to the nine SNRs currently available of the class 
(W28, W41, W44, W49B, W51C, Cygnus Loop, IC443, CTB 
37A, and G349.7+0.2) and reproduce the GeV- TeV spectral 
patterns as aforesaid. 

The basic parameters of each SNR studied here are 
listed in Table f with -E S nr = f0 51 erg and r\ = f0% as- 
sumed. Since they are all interacting SNRs, we only consider 
the MC which is in touch with the shock surface (i.e., the 
'contact' case, although in some cases not all the 7-rays is 
coming from a contact MC). Thus we adopt the value of the 
lower-limit of the integration of Eq. @ just the same as the 
value of the current shock radius, i.e., R c — AR c /2 = R s . 

In the 7-ray spectral fit, the parameters are adjusted 



in the following way. Given a certain cloud mass (chiefly 
adopted from observation), the 7-ray flux can be roughly 
determined by changing the diffusion correction factor \ 
(e.g., for point-like injection case, F-, oc \ ~ 3 ^ 2 (Gabici et 
al. 20f0)). The broadness of the platform is fitted by adjust- 
ing the cloud thickness AR C for the fixed \ and the slope of 
the "platform" and the slope above the break together with 
the TeV data can in turn be obtained with proper values of 
p and 5. 

3.1 G6. 4-0.1 (W28) 

W28 is one of the prototype thermal composite (or mixed- 
morphology) SNRs, characterized by shell-like radio emis- 
sion and centre-brightened thermal X-rays. It is an evolved 
SNR in the radiative phase, with an age estimate ranging 
from 35 to f 50 kyr (Kaspi et al. f993). In 7-rays the strongest 
TeV source is located at the northeastern edge of the rem- 
nant and positionally coincident with the well studied MC 
with which the remnant shock interacts (Aharonian et al. 
2008a). In this 7-ray emitting region, there are some of the 
OH masers points detected in W28 (Frail et al. f994). Re- 
cently, Fermi LAT observation found that the GeV source 
f FGL Jf80f.3-2322c coincides with the TeV source and thus 
demonstrates a broken power-law for the 7-ray spectrum 
with a break at ~ f GeV and photon indices of ~ 2.09 below 
the break and ~ 2.74 above the break (Abdo et al. 2010c). 
Three molecular clumps (Aharonian et al. 2008a) located 
outside the southern boundary of the remnant have differ- 
ent GeV and TeV brightness. In Paper I, we have explained 
the 7-ray spectra of the four sources via hadronic process 
of the diffusive protons, assuming different distances for the 
four sources. Here we also fit the spectrum of the north- 
eastern source, taking the cloud volume into account. The 
result is very similar to that in Paper 1. Unlike the sources in 
other SNRs studied in this paper, the "platform" below the 
break is unapparent (Figure (2j a)), entailing a large value of 
X and a small cloud thickness AR C . The evaluation of \ a l so 
sensitively depends on the cloud mass. 

3.2 W41 (G23. 3-0.3) 

HESS J1834-087 used to be representative of a class of 
"dark" 7-ray sources, which was not thought to be firmly co- 
incident with any detected sources (Aharonian et al. 2006). 
New H.E.S.S. analysis has revealed that the TeV emission 
can be divided into two components: the compact and ex- 
tended ones. While the TeV compact component compatible 
with pulsar candidate does not seem to have GeV coun- 
terparts, the extended one, with index 2.7, may spatially 
match the GeV emission detected by Fermi-LAT from the 
SNR W41 region (Mehault et al. 2011), making this SNR 
a promising candidate for generating the observed extended 
7-rays. The 7-ray spectrum in Fermi regime can be fitted 
by a single power-law with index 2.1, which causes a break 
between GeV and TeV. 

This field is also observed in radio with VLA, CO, and 
X-rays with XMM-Newton. The high-resolution CO obser- 
vation reveals its association with a giant molecular cloud 
with mass > 1O 4 M (Tian et al. 2007). The strong ab- 
sorption in X-rays together with the high the 7-ray to X- 
ray luminosity ratio (~ 11) indicates that the 7-rays may 



4 H. Li et al. 



(a) 




— — AR=Op<s : 




AR=8pc 






AR=10pc 












X=0.01; p=2.2; 5=0.5 \ 



10 100 1000 10000 100000 

Proton Energy (GeV) 




10 100 1000 10000 100000 

Proton Energy (GeV) 



E 

o 

> 

LU 



(b) 






p=1.8: 

— p=2.o: 








P=22; 














X=0.01; 5=0.5 





10 100 1000 10000 100000 

Proton Energy (GeV) 




100 1000 10000 100000 

Proton Energy (GeV) 



Figure 1. Spectra of the diffusive protons in the MCs for Rsnk — 20pc, t = 27000 yr, n = 1cm -3 , -Esnr — 10 51 erg, and R c = 25pc. 
In panel (a), AR C = Opc (blue), 8pc (black), and lOpc (red) are adopted, respectively. In panels (b)-(d), the black curves represent the 
case for AR C = 8pc and the red curves for AR C = lOpc. 



Table 1. Dynamical evolution parameters for interacting SNRs 



SNR 


age(kyr) 


distance(kpc) 


angular size 


Ksnr(pc) 


Reference 


G6.4-0.1(W28) 


42 


1.9 


48' 


13.3 


1, 2, 3 


G23.3-0.3(W41) 


100 


4.2 


33' 


20.2 


4 


G34.7-0.4(W44) 


20 


2.8 


30' 


12.2 


5,6 


G43.3-0.2(W49B) 


2 


10 


3' 


4.4 


7 


G49.2-0.7(W51C) 


30 


6 


30' 


26.2 


8,9 


G74.0-8.5(Cygnus Loop) 


20 


0.54 


180' 


14.1 


10,11,12 


G189.1+3.0(IC443) 


20 


1.5 


45' 


9.8 


6 


G348.5+0.1(CTB 37A) 


30 


11.3 


15' 


24.7 


13,14 


G349.7+0.2 


2.8 


23 


2' 


6.7 


14,15 



References.-(l) Paper I; (2) Velazquez et al. 2002; (3) Green 2009; (4) Tian et al. 2007; (5) Abdo et al. 2010a; (6) Seta et al. 1998 (7) 
Brogan & Troland 2001; (8) Koo et al. 2005; (9) Moon & Koo 1997a; (10) Green 1990; (11) Blair et al. 2005; (12) Miyata et al. 1994; 
(13) Sezer et al. 2011; (14) Reynoso & Mangum 2000; (15) Slane et al. 2002 



be of hadronic origin from energetic protons colliding with 
shocked GMC (Tian et al. 2007; Yamazaki et al. 2006). 



The GeV "platform" of this source is the broadest (~ 
l-100GeV) among the SNRs studied here (see FigureHfb)). 
It thus requires a low diffusion coefficient (x J$ 10~ 2 ) and 
large thickness of cloud (AR C comparable to R s ). With the 
parameters used (see Table 2), the well-determined TcV 
slope is also reproduced. 



3.3 W44 (G34.7-0.4) 

W44 is also a paradigm of "mixed-morphology" SNR, with 
semi-symmetric shape in radio with ~ 30' in size (Rho & 
Petre(1998)). It is considered to interact with a dense molec- 
ular cloud, with evidences including OH maser (Hoffman et 
al. 2005), molecular line broadening, CO line ratio (Castel- 
letti et al. 2007), etc. 

The 7-ray spectral energy distribution of W44 repre- 
sents a broken power-law, which breaks at 1.9GeV, with 
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Table 2. Fitted parameters for the 7-rays of interacting SNRs 



SNR 


P 


5 


X 


A R(pc) 


M cl (1O 4 M ) 


G6.4-0.1(W28) 


2.2 


0.5 


0.05 


0.5 


5 a 


G23.3-0.3(W41) 


2.1 


0.6 


0.004 


17 


7" 


G34.7-0.4(W44) 


2.4 


0.7 


0.02 


7 


3 C 


G43.3-0.2(W49B) 


2.4 


0.6 


0.03 


3 


1'' 


G49.2-0.7(W51C) 


2.2 


0.5 


0.04 


7 


l e 


G74.0-8.5(Cygnus Loop) 


2.2 


0.7 


0.04 


8 


0.025-f 


G189.1+3.0(IC443) 


2.3 


0.65 


0.01 


7 


0.5 9 


G348.5+0.1(CTB 37A) 


2.1 


0.6 


0.02 


5 


5.8 h 


G349. 7+0.2 


2.1 


0.7 


0.02 


3 


1.2* 



References.-(a) Aharonian et al. 2008a (b) Seta et al. 1998; (c) Seta et al. 2004; (d) Assumed; (e) Koo & Moon 1997 a,b; (f) Assumed 
(g) Cornett et al. 1977; (h) Reynoso & Mangum 2000; (i) Reynoso & Mangum 2001 



photon indices 2.06 at low energy and 3.02 at high energy. 
Most of the 7-ray emission is inferred to arise from the SNR 
other than from the unresolved pulsar based on the resem- 
blance between the 7-ray and infrared ring morphologies 
(Abdo et al. 2010a). Most recently, AGILE observation has 
also been reported with the 7-ray spectral down to 50 MeV. 
(Giuliani et al. 2011) 

Figure [21c) shows the fitting results of W44 using pa- 
rameters listed in Table 2. Our fitting not only well repro- 
duce the GeV "platform" spanned from 0.3 to 1.9 GeV, but 
also the spectrum with energy as low as 50 MeV, which 
is a robust advantage for hadronic origin of the 7-rays. It 
is noteworthy that the photon index above the break en- 
ergy (3.02) is the largest among the 8 studied SNRs; we 
adopt the upper limit (0.7) of the usual range for 5 so as for 
p(=2.4) not to deviate too much from the common values 
(1.8-2.2). Although the VHE spectrum of W44 is the very 
steep, we predict that the flux around TeV energy is slightly 
higher than the H.E.S.S. detection limit and much higher 
than the expected sensitivity of new generation TeV tele- 
scope Cherenkov Telescope Array (CTA) (de Cea del Pozo 
et al. 2009). 

3.4 W49B (G43.3-0.2) 

W49B is a "mixed-morphology" SNR at a young age. Near- 
infrared shocked H2 emission suggests that W49B is an inter- 
acting SNR (Keohane et al. 2007). Particularly high abun- 
dances of hot Fe and Ni, and relatively metal-rich core and 
jet regions are interpreted as evidence that W49B is origi- 
nated inside a wind-blown bubble inside a dense molecular 
cloud (Keohane et al. 2007). 

Recently, Fermi LAT reported a GeV source spatially 
coincident with W49B (Abdo et al. 2010d). The spectrum, 
again, exhibits a broken power-law shape with break energy 
~ 4.8GeV and the simple power-law can be rejected at a 
significance of 4.4<r. More interestingly, the best-fit spatial 
position of the LAT source which seems coincident with the 
brightest part of the Spitzer IRAC 5.8(im image strongly 
supports the scenario that the GeV emission is linked with 
the shocked gas. Recently, the H.E.S.S. collaboration has 
reported its detection in TeV and the photon spectrum 
is smoothly connected with GeV data points (Brun et al. 
2010). 

Due to the lack of definite value of the molecular mass 



around this SNR, we assume M c \ ~ 10 4 Mq for fitting the 
7-ray spectrum. The spectral slope below the break energy 
(< 4.8GeV) is relatively large (~ 2.2) and, consequently, a 
large value of p (= 2.4) is inferred. 

3.5 W51C (G49.2-0.7) 

W51C is also a member of "mixed-morphology" type of 
SNRs. A partial shell ~ 30' in diameter with a breakout fea- 
ture in northern part appears in the radio continuum (Moon 
& Koo 1994) . A shock-MC interaction was shown by the ob- 
servations of the shocked atomic and molecular gases (Koo 
& Moon 1997 a,b). 

In the 7-rays, the H.E.S.S. collaboration has detected 
an extended VHE source coincident with W51C (Fiasson et 
al. 2009), and the Milagro observation reported a possible 
excess in this direction (Abdo et al. 2009). Recently, Fermi 
LAT detects the flux in this field and shows that the GeV 
spectrum has a "platform" below the break energy ~ 2GeV 
and cannot be described by a single power-law, which makes 
it difficult for leptonic models to reproduce the GeV spec- 
trum (Abdo et al. 2009). Fang & Zhang (2010) explain the 
GeV- TeV spectrum with hadronic process, using the non- 
linear acceleration by low Mach number shock, and in this 
scenario the 7-rays arise from the acceleration region. 

Considering the protons diffusion into the adjacent MC, 
we perfectly reproduce the GeV- TeV spectrum with typical 
parameters (e.g., p = 2.2, 8 = 0.5) as shown in Figure [He). 

3.6 Cygnus Loop (G74.0-8.5) 

The Cygnus Loop is one of the most famous and well-studied 
middle-aged SNRs. The 7-ray emission mapped by Fermi- 
LAT is along the SNR limb and appears to be broken into 
four parts, three of which are essentially in correspondence 
with Ha bright patches (Katagiri et al. 2011). Actually, there 
have been some evidence of the presence of dense gas along 
the SNR boundary. In the west, CO emission in —25 — h30 
km s _1 interval seems to be adjacent to the SNR (Dame 
2001, see Katagiri et al. 2011) and two CO emission clumps 
at ~ +9.2 km s _1 (with a mass ~ IOOMq) seems to be well 
in contact with the optical filament (Scoville et al. 1977). 
Shock excited near-IR H2 line emission associated with the 
SNR is detected from the northeastern boundary (Graham 
et al. 1991a, b). Also, according to the X-ray studies, dense 
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wall of cavity is suggested to be along the eastern edge (Lev- 
enson et al. 1997; Levenson, Graham & Snowden 1999) and 
shock-cloud interaction is thought to be taking place there 
(e.g., Miyata & Tsunemi 2001; Zhou et al. 2010). Such a 
molecular-rich environment may act as a reasonable site for 
the hadronic origin of 7-ray emission from shock accelerated 
CRs, which is pointed out by Katagiri et al. (2011). 

Fermi LAT observation in the energy band 0.2-100 GeV 
shows a spectral break in the range 2-3 GeV. The 7-ray 
luminosity is ~ 1 x 10 33 erg s _1 , which is lower than those 
of other GeV-emitting SNRs. The 7-ray spectrum can be 
reasonably fitted by our model. The fitting result is shown 
in Figure [2{f) with parameters listed in Table 2. To specify 
AR C , we adopt the inner and outer radii 0.7' and 1.6' of the 
7-ray ring-like structure (Katagiri et al. 2011). 

3.7 IC443 (G189. 1+3.0) 

IC 443 is one of the most thoroughly studied interacting 
SNRs, with rich evidence including the OH masers from the 
central and southeast regions (Hewitt et al. 2006), CO line 
ratio (Seta et al. 1998), molecular line broadening (Dickman 
et al. 1992), near-infrared shocked H2 emission (Rosado et 
al. 2007), etc. 

In 7-ray s, there are several sources coincident with IC 
443 as observed by EGRET (Hartman et al. 1999), AGILE 
(Tavani et al. 2010), and Fermi (Abdo et al. 2010b) in GeV 
and MAGIC (Albert et al. 2007) and VERITAS (Acciari 
et al. 2009) in TeV. The spectrum can be represented by a 
broken power law with slopes of 1.93 and 2.56 and with a 
break at 3.25GeV and show very broad "platform" below 
the break (Figure [^f)). 

The multi-band photon spectra, with only the EGRET 
data then available in GeV band, was explained by using a 
time-dependent model for non-thermal particles emitting in 
the acceleration region (Zhang & Fang 2008). Torres et al. 
used the point source diffusion model assuming two clouds 
(with volumes ignored) in the vicinity of IC443 to explain 
the GeV- TeV connection (Torres et al. 2010). 

As shown in Figure [2ff), the 7-ray spectrum of IC443 
is fitted assuming one cloud with geometric volume. The 
broad GeV "platform" is reproduced with a low diffusion 
coefficient (x ~ 10~ 2 ) and a large thickness of cloud re- 
quired. 

3.8 CTB 37A (G348. 5+0.1) 

SNR CTB 37A has a partial shell with a extended breakout 
to the south (Kassim et al. 1991). The 1720 MHz OH masers 
(Frail et al. 1996) and morphologically correspondent CO 
emission (Reynoso & Mangum 2000) demonstrate that this 
SNR is interacting with MCs. 

In high energy 7-rays , H.E.S.S. has detected TeV emis- 
sion in the field of CTB 37A. For the absence of X-ray syn- 
chrotron emission, it is unlikely that the leptonic scenario 
could explain these VHE emission and a hadronic origin 
would be favorable (Aharonian et al. 2008b). Recently, an 
unresolved source lies within the eastern shell of CTB 37A 
was detected by Fermi LAT, and the spectrum in Fermi 
regime, too, shows a "platform" below 4.2 GeV (Castro & 
Slane 2010). The GeV-TeV spectrum is well fitted (see Fig- 
ure[2jg)) with parameters listed in Table 2. 



3.9 G349.7+0.2 

G349.7+0.2 is one of the radio and X-ray brightest SNRs in 
the Galaxy, with an blowout morphology (Slane et al. 2002). 
It has been shown to be an interacting SNR, with numerous 
evidences including the OH masers (Frail et al. 1996), CO 
line ratio of different transitions, molecular line broadening 
(Dubner et al. 2004), near-infrared (IR) shocked H2 emission 
(Reach et al. 2006), etc. 

The Fermi LAT observation shows a bright 7-ray source 
in its field with a significance from the evaluation of the test 
statistic ~ 10cr. The spectrum exhibits a broken power law 
with break energy ~ 16GeV (Castro & Slane). As yet, no 
TeV excess has been reported in this direction. With the 
spectral fitting (see Figure (Hh)), we predict that the flux 
around TeV energy is slightly higher than the H.E.S.S. de- 
tection limit and much higher than the expected sensitivity 
of CTA. 



4 DISCUSSION AND SUMMARY 

In this paper, the accumulative diffusion model for CRs es- 
caping from SNR shock (Paper I) is refined by considering 
the finite-volume clouds in the vicinity of SNRs. The vari- 
ation of the proton spectra with different model parame- 
ters are also shown for exploring the parameter space. This 
refined model is applied to nine Galactic SNRs which are 
thought to be interacting with the ambient dense material, 
and their GeV spectral breaks/ "platform" s, together with 
available TeV data, are naturally fitted. 

For the four SNRs (W51C, W28, W44 and IC443), 
Ohira et al. (2011) obtain the momentum breaks of proton 
spectra and successfully explain the 7-ray spectral shapes. 
In our model, in addition to fitting the spectral shapes, the 
real differential fluxes in 7-rays are also calculated, with the 
cloud masses considered. For all nine SNRs, the parameter 
5, which gives the energy dependence of the diffusion coeffi- 
cient, is used in an interval from 0.5 to 0.7, which is in the 
commonly favored range 0.3-0.7. The difference of the spec- 
tral fit parameters between our model and Ohira et al.'s 
(2011) (5 = 0.19 — 0.62) may be primarily caused by two 
factors. The first factor is the different assumption of escap- 
ing spectrum of protons: in their model, the time evolution 
of maximum energy and the delta-function momentum of 
the escaping protons are assumed, and in our model, typi- 
cal power-law distribution for escaping protons is used (e.g., 
Aharonian & Atoyan 1996; Torres et al. 2008). The second 
factor is the relative positions between the shock and the 
7-ray emitting MC: in their calculation, the MC is placed 
on the wall of stellar wind bubble separated from the shock 
(namely, the 'non-contact' case), while in our calculation, 
the 'contact' case is considered. 

For fitting the GeV- TeV spectra of these interacting 
SNRs, the diffusion coefficients (see Table 2) are required 
to be much smaller (x ~ 10 -2 ) than the averaged Galac- 
tic value ~ 1 and are also different from the values used in 
Ohira et al. (2011), which are in a wide range from 0.01 to 1. 
Such small values of x have been noticed in observation by 
several authors (e.g., Fujita et al. 2009; Giuliani et al. 2010; 
Paper I; Gabici et al. 2010; Torres et al. 2011) when apply- 
ing the proton diffusion to the 7-ray sources in the vicinity 
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Figure 2. 7-ray spectral data of the eight interacting SNRs and the model spectra. The black square and red dot show the data points 
observed by Fermi and H.E.S.S., respectively. The sensitivity curves for MAGIC (green), H.E.S.S. (red), and CTA (black) are also shown 
in (c), (f) and (i). Panels from (a) to (h) represent the spectra fit of SNRs W28, W41, W44, W49B, W51C, Cygnus Loop, IC443, CTB 
37A, and G349.7+0.2, respectively. (See text in Sec. [3]} 



of SNRs. Theoretically, this phenomenon is also studied by 
Monte-Carlo simulation, in which the CR streaming insta- 
bility can generate Alfven waves. The waves strongly scat- 
ter the particles and make the particle diffusion in the ISM 
around the SNR remarkably slow, and the diffusion coef- 
ficient is therefore strongly suppressed (Fujita et al. 2010; 
Fujita et al. 2011). This kind of suppression is significant if 
the ambient ISM is well ionized for the cosmic rays around 
the SNRs, because the neutral damping of Alfven waves is 
not effective. 

For small x> the typical diffusion velocity for certain 
particle energy may be smaller than the shock velocity, and 
the low energy cosmic rays will be continuously overrun 
by the SNR and reaccelerated. When considering particle 
acceleration, one needs a small Bohm diffusion coefficient 
(-Dboiim) to ensure particles to be proximate to the shock and 
repeatly cross it efficiently. The reason why particles escape 
is that their diffusion coefficient in the ISM (D esc ) is much 
larger than the one near the shock, namely D esc 3> -DBohm- 
In our paper, even if the diffusion coefficient is about two 
orders of magnitude smaller than the Galactic average, this 
condition is satisfied well. Even if some of the escaping par- 
ticles will be overrun by the shock, these particles only fill a 
much smaller volume near the shock surface than the volume 



of the particles diffusion region. Therefore, the probability 
of reacceleration process is quite small. 

Our accumulative diffusion model proves successful in 
explaining the 7-rays from the nine interacting SNRs, which 
is again indicative of the contribution of the protons escaping 
from SNR shocks to the diffusive CRs. This model shows 
its advantage in treating the hadronic emission which arises 
from the clouds near the SNR shocks. As a matter of fact, if 
the distance between the SNR and MC is considerably large 
than the SNR size, our result will be similar to the point-like 
injection case. 
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